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ABSTRACT 

Context. Massive star-forming regions are characterized by intense ionizing fluxes, strong stellar winds and, occasionally, supernovae 
explosions, all of which have important effects on the surrounding media, on the star-formation process and on the evolution of young 
stars and their circumstellar disks. We present a multiband study of the massive young cluster NGC 661 1 and its parental cloud (the 
Eagle Nebula) with the aim of studying how OB stars affect the early stellar evolution and the formation of other stars. 
Aims. We search for evidence of: triggering of star formation by the massive stars inside NGC 6611 on a large spatial scale (~ 10 
parsec) and ongoing disk photoevaporation in NGC 661 1 and how its efficiency depends on the mass of the central stars. 
Methods. We assemble a multiband catalog of the Eagle Nebula with photometric data, ranging from B band to 8.0 yum, and X-ray 
data obtained with two new and one archival CHANDRA/ ACIS-I observation. We select the stars with disks from infrared photometry 
and disk-less ones from X-ray emission, which are associated both with NGC 6611 and the outer region of the Eagle Nebula. We 
study induced photoevaporation searching for the spatial variation of disk frequency for distinct stellar mass ranges. The triggering of 
star formation by OB stars has been investigated by deriving the history of star formation across the nebula. 

Results. We find evidence of sequential star formation in the Eagle Nebula going from the southeast (2.6 Myears) to the northwest 
(0.3 Myears), with the median age of NGC 661 1 members ~1 Myear. In NGC 661 1, we observe a drop of the disk frequency close 
to massive stars (up to an average distance of 1 parsec), without observable effects at larger distances. Furthermore, disks are more 
frequent around low-mass stars (< 1M Q ) than around high-mass stars, regardless of from the distance from OB stars. 
Conclusions. The star formation chronology we find in the Eagle Nebula does not support the hypothesis of a large-scale process 
triggered by OB stars in NGC 6611. Instead, we speculate that it was triggered by the encounter (about 3 Myears ago) with a giant 
molecular shell created by supernovae explosions about 6 Myears ago. We find evidence of disk photoevaporation close to OB stars, 
where disks are heated by incident extreme ultraviolet (EUV) radiation. No effects are observed at large distances from OB stars, 
where photoevaporation is induced by the far ultraviolet (FUV) radiation, and long timescales are usually required to completely 
dissipate the disks. 



1. Introduction 

In the first decade of the 21th century, a large number of ac- 
tive star-forming sites in our Galaxy have been identified and 
studied in detail with a multi-wavelength approach, allowing to 
advance our understanding of star formation (see, for example, 
Feig elson & Montmerld 119991) . even if we yet lack a full un- 
derstanding of how molecular clouds turn into stellar clusters. 
In particular, it is still unclear what initiates the cloud collapse 
and what the environmental effects are on the stellar population 
formed after the collapse. 

The initial impulse for cloud collapse may be provided by 
shocks created by intense stellar winds from a large amoun t 
of massive stars or supernovae explosions dCho & KangL 12007). 
Sometimes they can form giant molecular shells, which can be 
associated to star-formation events, such as the Barnard loop and 
the E ridanus loop in the Orion complex dRevnolds & Q gden, 
[1971 . 

It is a debated question whether the stellar population pro- 
duced by the collapse of a molecular cloud depends on the en- 
vironment in which stars form. This question is also related 
to the universal relation between the cluster age and the frac- 
tion of members that still bear a disk (i.e. the disk frequency, 



Send offprint requests to: mguarce@astropa.unipa.it 



lHaisch et all l200li) . In particular, the role played by OB stars 
in young clusters is unclear. Their energetic radiation, which 
ionizes the HI clouds into HII, and their intense winds disperse 
the cloud, which halts the star-formation process. On the other 
hand, the high pressure provided by their stellar wind and en- 
ergetic radiation may compress the clouds, inducing their ra- 
diative im plosion and the formatio n of subsequent generations 
of stars dMcKee & Ostriken, 120071) . The intense UV radiation 
from massive stars may also dissipate the protost ellar cores 
emerg ed from the evaporating parental clouds dHester et all 
1996) and the nearby circumste llar disks by photoevaporation 
(i.e. lStorzer & Hollenbachlll999t) . which affects the mass accre- 
tion onto protostars and the formation of planetary systems. 

The effects of massive stars on the evolution of Young Stellar 
Objects (YSOs) and on the star-formation process can be ana- 
lyzed with several methods. For instance, for young, not yet re- 
laxed clusters, the study of the sequential star formation from the 
positions of massive stars outward can give a clue about the rel- 
evance of induced cloud-compression on the formation of new 
generations of stars. This approach has been succ essfully applied 
in sev eral cases, such as the HII regions RCW 7 9 dZavagno et all 
2006), Sh 217, Sh 219 dDeharveng et alll2003l). as well as i n the 
star-formin g regions 30 Dora dus dWalborn & Bladesl l 1997b and 
NGC 1 893 dNegueruela et all 120071 7" 



1 



M. G. Guarcello et al.: Disks frequency and age in M16 



The outcomes of the investigation on disk frequencies and 
Initial Mass Function (IMF) in young clusters are much more 
debated. It is well known that the d isk frequency decre ases on 
average with increasing cluster age (Haisch et all 120011) . There 
is no general consensus yet whether the observed disk frequen- 
cies are ind ependent of the cl uster population of OB stars (as in 
NGC 661 l. lOliveira et al.l2005h or if instead they decrease close 
to OB stars as a conseque nce of externally i nduced photoevap- 
oration (as in NGC 2244. iBalog et al.ll2007h . Furthermore, the 
universality of the IMF seems to imply that the dissipation of 
protostellar cores by the UV radiation from massive stars does 
not significantly affect the mass accretion onto YSOs. 

We investigate these open issues for the Eagle Nebula (Ml 6) 
and the associated young cluster NGC 6611. In Sect. [2] we re- 
view some of the main characteristics of both the Eagle Nebula 
and NGC 66 1 1 , mostly focusing on the morphology of the re- 
gion and on the existing evidence about the feedback by massive 
stars on the star formation and disk evolution. Since a multi- 
wavelengths approach is necessary to select and classify T-Tauri 
stars and to obtain information about their stellar parameters, in 
this work we analyze a set of optical, infrared, and X-rays data, 
which are described in Sect. [3] These data are combined in in 
a multi-band catalog of the region, which is described in Sect. 
|4] Appendix [C] is devoted to the procedure used to match all the 
data and obtain the multi-band catalog. The selection of stars 
with infrared excesses due to a circumstellar disk, performed 
with suitable reddening-free color indices, optical and infrared 
color-magnitude diagrams and the IRAC color-color diagram is 
described in Sect. [5] In that section we also analyze the contam- 
ination of the list of selected candidate stars associated to M16 
by both foreground and background objects. The final member 
list is summarized in Sect. [6] The effects of the UV radiation 
that is emitted by massive stars on the evolution of circumstel- 
lar disks in M16 are analyzed in Sect. [7] where we also com- 
pare the timescales for externally induced disk photoevaporation 
and "normal" dissipation at different masses of the central star. 
The triggering of star formation in M16 by the massive stars 
in NGC 661 1 is analyzed in Sect. [8j where we find finding the 
chronology of star formation across the cloud. All the results are 
discussed in detail in Sect. [9] In Appendix [A] we describe the 
stellar content of three peculiar regions of the Eagle Nebula; in 
Appendix [B] we derive the extinction map of the cloud; while 
in Appendix [C] we describe the adopted catalog cross-matching 
procedure. 

2. The Eagle Nebula (M16) and NGC 6611 

The Eagle Nebula, one of the most active star-forming sites in 
the Sagittarius arm, and the young open cluster NGC 6611 in 
its center are ideal targets to study star formation in the pres- 
ence of massive stars. The distance to the cloud has been eval- 
ua ted by several authors ; here, we adopt the value computed 
in iGuarcello et all d2007l) . hereafter GPM07, which is equal to 
1750parsec. 

Starting from the study of Ide Winter et al ] (119971) . several 
works have been devoted to the massive stars of NGC 66 1 1 . The 
cluster contains a remarka ble population of OB stars (54, with 13 
O stars, Hillen brand et all 1 19931) . which are concentrated mostly 
in the central cavity of the molecular cloud, which is about 5.5' 
(2.2 parsec) wide. These massive stars are responsible for an in- 
tense UV field in the cluster center; for comparison, the three 
most massive stars in NGC 661 1 alone emit about five times the 
ionizing flux of the O stars in the center of the Trapezium in 
the Orion Nebula Cluster (ONC), where disk photoevaporation 



has been directly observed (IStorzer & Hollenbachi Il999t) . The 
most massive star in the cluster is W205 (a =18*:18 m :36*.48, 
6 = -13 rf : 48 m :02'.5) with a mass of 75-80M o and a sp ectral class 
03-05V (iHillenbrand et allll993blEvans et alll2005l) . The main 
source of ionizing fl ux in the nebula is W 205, accounting for half 
of the total UV flux dHester et all[l996h . 

Only few OB stars are far away from the cluster center. 
The most noticeable among them is W584, an 09V star at 
a =18:18:23.69, 6 =-13:36:28.1, about 6.5' from the approxi- 
mate cluster center in the northwest direct ion. As inferred by the 
[8.0] IRAC image of the Eagle Nebula (IGuarcello et all 12009. 
hereafter GMD09), this star should have created a cavity in the 
nebula of about 1.5' of radius (0.7 parsec). 

Several studies of the Eagle Nebula concern the trigger- 
ing of star formation by the massive members and their influ- 
ence on the evolution o f both circumstellar disks and YSOs. 
IHillenbrand et alJ (1 19931) found that the massive population of 
NGC 6611 has a mean age of 2 Myears, with a spread of 1 
Myear, even if one of the massive members (with a mass of 
30M o ) is about 6 Myear old. Subsequent studies in the central re- 
gion of NGC 661 1 discovered a large age spread also for the low- 
mass members, w hose age ranges from les s than 1 Myear to ~ 3 
Myears (GPM09, Indeb etouw et al.ll2007h . The reasons for this 
age spread and the sequenc e of sta r formation are still unclear, 
but iGvaram adze & Bomans (2008) hypothesized that the oldest 
massive stars are blue stragglers, while the younger pre-main se- 
quence (PMS) stars are part of distinct generations formed in 
successive star-formation episodes. It is reasonable to assume 
that the oldest massive stars triggered the star formation in the 
central region of NG C 6 61 1 in the last 2 Myears. However, both 
iBeUkovet aT] d2000t) and llndebetouw et all d2007l) found no clue 
for triggered star-formation in the central region of the cluster. 

Even if no proof of sequential star formation has been found 
in NGC 661 1, there are indications of a short spatial scale feed- 
back of OB stars on the star formation and the evolution of 
the cloud. For instance, the radiation from OB stars continu- 
ously models some cloud structures, such as the well known 
"Pillars of Cr eation", which con t ain deeply embedded disk- 
bearing stars (| H ester et all 1 1996b iMcCaughrean & Andersenj 
120021: iHealv et all 12004 ISugitani et alll2007h . 

Other active star-forming sites are present in the Eagle 
Nebula. The Bright Rimmed Cloud SFO30 is a region bright 
in IRAC bands, ~ 6' away from the cluster center (~ 3 parsec) 
northward, where water maser and sour ces with excesses in near- 
infrared (NIR) bands have been found (|Hillenbrand et"ai1 11993b 
IHealv et alj|2004t llndebetouw et alj|2007l GMD09). An embed- 
ded massive star in for mation has been found 8 ' (3.8 parsec) 
west of the cluster center (llndebetouw et alll2007h . In the oppo- 
site direction (~ 7' eastward from the cluster center, 3.3 parsec) 
there is another elongated structure called Column V, probably 
modeled by the radiation from the cluster massive stars, where 
star-f ormation activity is proved by water masers (|Healv et all 
2004) and can didate Herbig-Haro ob jects (Meaburn & Walsh, 
1986). Finally, Indebetouw efal] d2007l) found a new cluster of 
embedded objects inside a de nse region of the Eagle N ebula, 16' 
northeast of the cluster center. IGuarcello et a l. (2009) found that 
due to the strong absorption, only one of these young stars has 
a faint optical counterpart. Some of these active star-formation 
sites will be discussed in more detail in AppendixlAl 

The effects of massive stars on star formation and disk 
evolution in NGC 6 611 have been studied by other authors: 
lOliveira et alJ (l2005h evaluated the disk frequency in the cen- 
tral region of the cluster to be equal to 58%, comparable with 
that found in clusters with smaller populations of massive stars; 
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lOliveira et "all (120091) obtained the IMF of candidate cluster 
members, which has the same slope of the IMF as that of the 
Taurus region. These studies suggest that in NGC 661 1 massive 
stars have no evident effect on the evolution of the other cluster 
members. On the other hand, GPM07 and GMD09 found that 
the disk frequency in the central region of the cluster decreases 
at small distance from massive stars, with a minimum value of 
~ 16% in the regions where the disks are illuminated by intense 
UV flux. Their results are compatible with a scenario in which 
UV radiation from massive stars destroyed the nearby circum- 
stellar disks by induced photoevaporation. 

iGuarcello et alj20"09l compiled a list of 1264 stars associated 
with M16: 474 ClassI/II objects selected for their infrared ex- 
cesses in 2MASS and IRAC bands in the 33' x 34' field centered 
on NGC 661 1 and 790 ClassIII PMS stars selected for their X- 
ray emission in the 17' x 17' region centered on NGC 6611. In 
this paper we extend the selection of disk-less members to the 
entire Eagle Nebula region. Catalogs of the whole region have 
been published by Hillenbrand et al 
BVIJHK photometry), Belikov et al 



1993 



1999 



(1026 sources with 
(2200 sources with 
BVJHK photometry) and GPM07 (38995 with BVIJHK and X- 
ray data). 




□ashed boxes: WFCAM fields 

Fig. 1. WFI-I image of the Eagle Nebula in / band (the 4x2 
mosaic). The fields covered by UKIRT/WFCAM (dashed boxes) 
and ACIS-I (inclined boxes) are also shown. NGC 6611 lies in 
the center of the ACIS-I C - field. North is up, East is on the 
left. 



3. Data analysis 

We compiled an optical-infrared-X-rays catalog of sources in the 
33' x 34' region centered on NGC 66 1 1 . The main improvements 
with respect to the catalog published by GPM07 are the use of 
new X-ray data in the outer regions of the nebula and the exten- 
sion to fainter magnitudes in the infrared JHK bands, taking ad- 
vantage of the UKIDSS/GPS survey (see Sect. [372t . Moreover, 
the instrument used in the UKIDSS surveys has an improved 
spatial resolution compared to that of the 2MASS survey. This 
is crucial in studies of crowded regions such as NGC661 1. 

In summary, optical data in BVI bands have been obtained on 
2000 July 29 with WFI/2.2m@ESO as part of the ESO Imaging 
Survey (EIS, iMomanv et al.L 1200 lb . and they have been re- 
duced and calibrated in GPM09. Infrared data in th e JHK bands 
are ta ken from the Galactic Plane Survey (GPS, iLucas et ail 
2008) of the United Kingdom Infrared Deep Sky Survey 
(UKIDSS) , based on observa tions with the Wide Field Camera 
(WFCAM. ICasali et alll2007h on the United Kingdom InfraRed 
Telescope (UKIRT). Spitzer/IRAC data are taken from the cat- 
alog Galacti c Legacy Infrared Mid -Plane Survey Extraordinaire 
(G LIMPS E, Benja min et all l2003). and they have been analyzed 
in llndebetouw et all d2007l) and GMD09. Both UKIDSS and 
IRAC data cover the complete WFI field of view (FoV), with 
the exception of a small region in the northwest for the IRAC 
data. X-ray data of the cente r of the Eagle Nebula are taken 
from the CHAN DRA/ACIS-I (IWeisskopf et all 120021) observa- 
tion presented in Linsk v et al .1 (120071) : in addition, we present 
here two new CHANDRA/ACIS-I observations pointed toward 
east and northeast with respect to the cluster (PI: Guarcello). 
Table Q] resumes the three CHANDRA observations. The cen- 
ter of the new observations have been chosen in order to have 
ColumnV and the embedded cluster (Sect. [2} near the centers 
of the ACIS FoV, where the PSF is narrow and the sensitivity 
is better. Thereafter, we will call these fields C-field, E-field and 
NE-field, respectively. 

Figure [T] shows the 33' x 34' image in / band of the Eagle 
Nebula, obtained with WFI. The FoVs of WFCAM and ACIS- 
I are also shown, marked with the dashed and inclined boxes 
respectively. 



3.1. X-ray data 

The reduction of X-ray data was performed with CIAO 4.03 and 
the CALDB 4.1.1 calibration files. In order to perform a consis- 
tent study, we analyzed both the new observations in the outer 
fields and that in th e C-field, which was already analyzed in 
iLinskv et all d2007l) . For each observation, a Level 2 event file 
was produced with the ACIS .PROCESS -EVENT CIAO tool from 
the Level 1 event file. In the Level 2 files we retained only events 
which produced a charge distribution in the detector that is in- 
compatible with cosmic rays (grades -0,2, 3, 4, 6 and status-0). 
The 0.5" events position randomization added in the Chandra 
X-ray standard data processing was removed, because it could 
have affected the sources position determination. 

The 1 keV exposure map, necessary for point-source detec- 
tion, was produced with the MERGE J^LL task. Point-source de- 
tection was performed with PWDetecffi, a wavelet-based detec- 
tion algor ithm developed at IN AF-Osservatorio Astronomico di 
Palermo dDamiani et all H997I) . PWDetect found 1158 sources 
in the C-field, 363 in the E-field and 315 in the NE-field, with 
a threshold limit of 4.6<x, corresponding to 10 expected spuri- 
ous detections, and considering only photons in the 0.5-8 keV 
energy band. This 4.6cr value has been adopted because the ex- 
pected amount of spurious detections is a negligible percentage 
with respect to the sources total number, and it is smaller than the 
number of sources we would lose with higher <x values. After a 
careful visual inspection of the detected sources, we removed 
1 1 double detections of the same source found with wavelets at 
different spatial scales. Using WEBPIMMSQ and the relation Lx 
vs. stellar mass found by Preibis ch et all d2005l) . we evaluated 
the detection limit in the C-field to be just 0. 1 solar mass deeper 
than the corresponding limit in the external ACIS-I fields (from 
0.3 M e to 0.2 M Q ). A detailed account of the properties of X-ray 
sources we detected will be the subject of a forthcoming paper. 
The full list of X-ray sources, compiled with the procedure ex- 



1 http://cxc/harvard.edu/ciao 

2 http://www.astropa.unipa.it/progetti_ricerca/PWDetect 

3 http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 
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Table 1. CHANDRA/ AGS -I observations of the Eagle Nebula. The E-field has been observed on four different days. 



Obs. ID 


Field name 




RA 
(J2000) 


DEC 
(J2000) 


Exposure time 
Ksec 


Date 




978 


C- field 


18 


18 : 44.72 


-13 : 47 : 56.50 


78 


2001 - 07 


-30 


8931 


NE-field 


18 


19 : 12.00 


-13 : 33 : 00.00 


80 


2008 - 05 


-28 


8932 


E- field 


18 


19 : 36.00 


-13 : 47 : 24.00 


31 


2008 - 06 


-02 


9864 


E- field 


18 


19 : 36.00 


-13 : 47 : 24.00 


23 


2008 - 06 


-07 


9865 


E-field 


18 


19 : 36.00 


-13 : 47 : 24.00 


17 


2008 - 06 


-04 


9872 


E-field 


18 


19 : 36.00 


-13 : 47 : 24.00 


9 


2008 - 06 


-09 



plained in Appendix consists of 1755 X-ray sources in the 
three fields. 

3.2. UKIDSS data 

Near infrared data in JHK bands were extracted from the 
UKIDSS/GPS catalog. The GPS survey covers the entire north- 
ern and equatorial Galactic plane ar ea of the sky plus th e Taurus- 
Auriga-Perseus molecular complex dLucas et alll2008l) . and also 
includes the Eagle Nebula. The observations were made with the 
UKIRT/WFCAM telescope and consist of short exposures for a 
total integration time of 80 sec. (J and H bands) and 40 sec. (K 
band), respectively. 

The WFCAM is equipped with a 2 x 2 array of CCDs, 
2048 x 2048 pixels each, with a spatial scale of 0.4" per pixel, 
smaller than the 2MASS scale ( 1 " per pixel ); the typical FWHM 
is 0.8". Infrared observations with such a high spatial resolution 
are crucial for studies of crowded clusters such as the core of 
NGC 6611. This high spatial resolution reduces the possibility 
of confusion in merging the detections in different bands, which 
is important because spurious matches in this merging process 
can be mistaken as real stars with non-photospheric colors. 

We selected the UKIDSS sources in the WFI FoV with 
priorsec-0 (sources not duplicated in observations of adjacent 
fields), merged j:lass--l (star-shaped sources), JHKerrbit<256 
(sources with moderate quality warningsQ), for a total of 159999 
sources down to / = 19'". Note that 10% of sources with a 
high possibility of being noisy (merged jzlass-G) can be stars 
and that sources classified as galaxies or probable galaxies 
(merged jzlass-+l ,-3) can be unresolved pairs of point sources. 
These sources have been excluded from the catalog to produce a 
more reliable list of cluster members. 

We matched UKIDSS and 2MASS catalogs of the Eagle 
Nebula with a 0.3" matching radius (see Appendix [C]), find- 
ing 17396 sources in common between the two catalogs, with 
few multiple matches. The unmatched sources are 8522 in the 
2MASS catalog (mostly bright sources) and 142603 in the 
UKIDSS catalog (mostly very faint sources). The large number 
of stars detected only in the UKIDSS reflects the depth of this 
catalog. 

The zero-point calibration of t he GPS catalog wa s performed 
using 2MASS standard sources dLucas etal.L |2008). Still, using 
stars with good photometry in common between the 2MASS and 
the UKIDSS catalogs, we found small offsets between 2MASS 
and UKIDSS photometry: +0.06'" (J), -0.04 m (H), -0.02 m (K). 
Because magnitudes of UKIDSS/GPS were calibrated in the 
2MASS photometric system, we corrected the magnitudes of our 
UKIDSS sources for these offsets. We also evaluated the satura- 



4 http://surveys.roe.ac.uk/wsa/index.html 
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Fig. 2. J 2 mass -J ukidss vs. J ukwss for the stars with good pho- 
tometry. The horizontal line shows the small systematic offset; 
the dotted line is the adopted J saturation limit. 

tion limits of our UKIDSS catalog as the magnitude at which 
the differences between UKIDSS and 2MASS magnitude stars 
to deviate from the median. The case of the J band is shown in 
Fig-El where the horizontal line marks the median of the magni- 
tud e differences . The li mits are J = 12 m , H = 12.5'", K = ll m . 

iLucas et al .] (120081) claim that the millimagnitude errors 
quoted in the UKIDSS/GPS catalog are not reliable, and that 
they should be replaced with a lower limit of 0.0 2 1 ". In this work 
we follow the suggestion of ILucas et al . (2008), so that we do 
not underestimate the systematic uncertainties. 

4. Multiband catalog 

In the 33' x 34' region centered on NGC 6611, we selected 
28827 optical sources (down to V = 23 m ) 159999 UKIDSS 
sources (down to J = 19 m ), 41985 IRAC sources (down to 
[3.6] = 13 m ) and 1755 X-ray sources (down to F x = 1.5 x 10~ 15 
erg x cmT 2 x s~ l ). We referred the astrometry of these catalogs 
to that of 2MASS PSC, and we cross-correlated the sources in 
these catalogs with the procedure described in Appendix ICl 

Figure|3]shows the J versus J — K diagram for all stars in the 
WFI FoV with both crj < 0.1'" and ctj-k < 0.15 m . Hereafter, in 
all the color-color and color-magnitude diagrams shown in this 
work we will plot only objects which satisfy this color and mag- 
nitude selection. The region betwee n the 0. 1 and the 3 Myears 
isochrones (from ISiess et all |2000) defines the cluster locus, 
which is well traced by X-ray sources (marked with boxes). The 
region on the left of the cluster locus is populated by field MS 
stars. This locus is populated down to the sensitivity limit by 
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o.i u. 



1 2 3 4 5 6 7 

J-K 

Fig. 3. J versus J-K diagram for all stars in WFI FoV (points). 
The dashed lines are the 0.1, 0.25, 1, 3, 5, 10 Myears isochrones 
from ISiess et aTl d2000t) at the cluster distance and reddening 
(GPM 07); the reddening vector w as obtained from the reddening 
law of iRieke & Lebofskvl £[985). The mass values listed on the 
right side of the diagrams correspond to stars along the 1 Myear 
isochrone. The overplotted squares mark the X-ray sources. 

sources with both UKIDSS and WFI detection, while those de- 
tected also with IRAC (not shown in the figure) populate the 
field star locus only down to J ~ 15'", due to the shallowness of 
this sample of stars. The region on the right of the cluster locus 
is that of background giants. Few optical sources fall in this lo- 
cus, which is dominated by only UKIDSS sources down to the 
sensitivity limit. 

5. Candidate members with circumstellar disk 

5.1. YSOs from IRAC colors 

Candidate young stars with infrared excesses have been se- 
lected using two criteria: IRAC color-color diagram and suitable 
reddening-free color indices (Q indices). 

The selection of disk-bearing YSOs with the IRAC color- 
color diagram has been already presented in GMD09, who found 
143 ClassII stars, 13 ClassI, and 22 objects not clearly classified 
(ClassI or highly embedded ClassII stars). Figure 3 of GMD09 
shows the IRAC color-color diagram and the selection of YSOs 
associated with M16. 

Stars with excesses in the IRAC bands were also selected 
using the Q color indices, which are independent from inter- 
stellar extinction. Q indices have been described in detail in 
iDamiani et all d2006l) , GPM07 and GMD09. In particular, the 
Q indices used to select stars with excesses in IRAC bands, 
Qvu[sp], are defined as 

Qvmsp] = (V -/)-(/ - [sp]) x Ev-dEj. m , (1) 

where [sp] is the magnitude in an IRAC band and V — I 
is representative of photosp heric emission. In Eq. [D the red- 
dening Ey-i is obtained from lMunari & Carrarol (Il996l) . Ej-[ sp ] 



was inferred from the reddening law of Mathisj (Il990h . In the 
Qvu[sp] vs - (J - [sp]) diagram (see GMD09), stars with ex- 
cesses in the [sp] band have smaller indices than those of nor- 
mal stars, while the extinction vector is horizontal. This allows 
the separation of stars with excesses from extincted sources. In 
the Qvuyp] vs - (J ~ [ s p]) diagrams, stars with excesses are de- 
fined by us as those with the Q index smaller (below 3 x <tq) 
than the limit for the photospheric emission. We select 343 can- 
didate stars with disks considering Qvuisp) indices, which were 
mostly already found with the 2MASS photometry (the selection 
of stars with disks with UKIDSS and 2MASS will be compared 
in Sect. 15. 3b . The selection of stars with excesses in IRAC bands 
that we performed here with Q indices is only slightly differ- 
ent from that in GMD09, because we use here the UKIDSS J 
band, instead of the 2MASS J. Both selections share the same 
lower limit in J (~ 16'"), due to the completeness limit of IRAC 
data, but have a different upper limit (J ~ 13.5"' with Jukidss , 
J ~ 10 ffl with J 2 mass)- In GMD09 the selection of YSOs with 
the IRAC color-color dia grams is compared in detail with that of 
llndebetouw et al.1 d2007l) . 

The nature of the Qirac - excess stars is shown in the di- 
agrams in Fig. |4] In the J vs J - K diagram (left panel), the 
Qirac - excess stars lie in the cluster locus with J-K redder 
than the colors corresponding to the 0. 1 Myear isochrone, likely 
due to a K excess. These stars lie in the cluster locus also in the 
optical color-magnitude diagrams. Because Qvij[sp] requires the 
combination of WFI, UKIDSS and IRAC photometry, this selec- 
tion cannot be extended to the faint UKIDSS sources. The IRAC 
color-color diagram (right panel) shows that all these stars are 
ClassII objects, mostly with low reddening. 

5.2. YSOs from UKIDSS colors 

Stars with excesses in the UKIDSS bands have been selected 
with the Q indices defined in GPM09. Four indices compare V-I 
with NIR colors: 

Qukidss = (V-I)-(A-B)x Ev-i/Ea-b, (2) 

where A -B assumes the values I- J, J-H , J-K or H-K and 
the red deni ng Ea-b is calculated acco rding to lRieke & Lebofsky 
(Il985h and iMunari & Carrar o ( 1996). Another index compares 
J-H with H-K: 

Qjhhk = (J-H)-(H-K)x Ej- H /E H - K , (3) 

Because of detection in the optical bands is not required, this 
index allows us to select more absorbed stars than the Qvij[ S p] 
and the other Qukidss indices. Figure [5] shows the Qjhhk vs. 
H-K diagram for stars in the WFI FoV (gray points), which 
has been used to select stars with excesses in K band (marked 
with circles). These stars are those with a more negative index 
than the corresponding photospheric limit by more than 3<x G . 
The other Qvu[sp] an d Qukidss diagrams are analogous. 

Table [2] compares the number of candidate stars with ex- 
cesses in JHK bands selected with UKIDSS (this work) and 
2MASS photometry (GPM09); the number of selected stars is 
126^1 and 345, respectively. This large difference cannot be ex- 
plained only by the depth of the UKIDSS catalog. 



5 The number of stars selected with Qukidss is weakly depending 
on the adopted saturation limits (see Sec. 13. 2i ): increasing the limits 
by one magnitude we selected two more stars, decreasing them by one 
magnitude we selected 29 stars less. 
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Fig. 4. Color-magnitude diagrams for the stars (gray points) in the WFI Fo V ( J vs J - K in the left and the IRAC color-color 
diagram in the right panel). Squares mark the subsample of stars with excesses in IRAC bands selected with Qvij[sp] indices with 
good UKIDSS photom etry (left panel) and good IRAC photometry (right panel). In the left panel the ZAMS and isochrones are 
from lSiess et al.l J2000). The mass values listed on the right side of the diagrams correspond t o stars along the 1 Myear isochrone. 
In the right panel the box delimits the approximate locus of ClassII YSOs ([Allen et al.ll2004j) . The dashed lines in the right panel 
roughly separate reddened photospheres ( left), reddened ClassII YS Os (center), and Class I YSOs. The reddening vectors in the right 
panel with A v =30 m and A K =5 m are from lMegeath et all d2004 and lFlahertv et all (120071) . respectively. 



Table 2. Number of stars with excesses in JHK bands selected 
with the Q indices defined with UKIDSS (second column) and 
2MASS (third column) photometric data. 



Q index 


UKIDSS detection 


2MASS detection 


QjHHK 


448 


22 


QvilJ 


39 


34 


QviJH 


375 


91 


Q.VIJK 


751 


159 


Q.VIHK 


609 


136 



Figure|6]shows the V vs. V-I and J vs. J—K diagrams for all 
the stars in the WFI FoV (gray points) and the stars with excesses 
selected with the Qukidss indices. In the right panel a signifi- 
cant number of these stars are faint and highly reddened. In the 
left panel it is evident that a great number of Qukidss ~ excess 
stars lie in the locus of field stars instead of in the cluster lo- 
cus. It is possible that a physical mechanism, such as gas ac- 
cretion from the disk to the central star or scattering of stellar 
light into the line of sight by disk's grains, could produce these 
blue optical colors in young disk-bearing stars, which fall in 
the field locus in the various optical color-magnitude diagrams 
(iHartmann & Kenvon|[l990t iHillenbrandll 19971: iGuarcello et al l 
l2010h . These effects, however, affect only the optical colors, 
while most of these stars fall in the field star locus also in di- 
agrams defined with NIR photometry. This suggests that some 
stars have been wrongly selected with Qukidss indices. We 
decided, then, to confirm as disk-bearing cluster members the 
Qukidss ~ excess stars which share at least one of the following 
features: 

1. they are compatible with PMS stars in the IRAC color- 
magnitude diagrams (defined by two or three IRAC bands), 



i.on 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

. Norma^ Stars 
s 4 1 

J \ \ Recdened Stars 




Stars with excess o 8 o ° 

o o a Q 

-2.0 L i i i I , , , o_J , , , , ° ° , i i i i I 

0.0 0.5 1.0 1.5 2.0 2.5 

H-K 

Fig. 5. Qjhhk vs. H - K diagram for the stars in the WFI FoV 
(gray points). Dashed and dotted lines separate the loci of normal 
stars, reddened stars, and stars with excess (circles). The solid 
lines are the isochrones as in Fig. [3] 



in the V vs. V - I, I vs. / - J or J vs. J - K diagrams (424 
stars); 

2. they are also Qirac ~ excess stars (102 stars); 

3. they fall in the ClassII locus in the IRAC color-color diagram 
(49 stars). 

Adopting these criteria, we selected 575 candidate disk- 
bearing stars, down to 0.2 M & . 
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We discarded 399 Qukidss ~ excess stars from the list of 
candidate members because 

1. they he in the field-stars loci both in the optical and in the 
infrared color-magnitude diagrams (389 stars); 

2. or they have normal colors in the Qvij[sp] vs. {J - [sp]) dia- 
grams (10 stars); 

We have no firm explanation for the 399 objects for which 
the Q selection seems to fail. One possibility is that these stars 
are false matches between the optical and infrared catalogs. 
However, t he expected amount o f spurious identifications, eval- 
uated as in lDamiani et alj d2006l) . is 107, 3.7 times smaller than 
the amount of false Q selections. More likely, these selections 
are induced because the adopted reddening laws is not reliable, 
modifying the ratio of color excesses used in Eqs.[2]and[3] This 
should happen both because there can be problems in the calibra- 
tion of UKIDSS/GPS data on the 2 MASS photometric system, 
especially for very reddened objects dLucas et aUl2008h . and be- 
cause the adopted extinction laws are defined for the usual inter- 
stellar Ry = 3.1, while in the direction of NGC 6611 the red- 
dening law is anomalous (the value Ry = 3.1 is used for the 
ratio between color excesses in Q definitions, even if a slightly 
higher value, 3.27, has been found in GPM07). Different values 
of Ry change both the slope of the reddening vectors in the color- 
magnitude diagrams and the relations between the reddening in 
various colors. The reddening laws we used do not give the exact 
relations between the reddening at different bands if Ry differs 
from the usual 3.1. 

The remaining 290 stars are retained in the final catalog as 
"probable members", because we could not assess their nature, 
with a flag indicating their "uncertain members status". Most of 
them are faint QjHHK-excess sources without any other counter- 
part. We retained these stars because most of the new detections 
of low-mass disk-bearing stars provided by UKIDSS should be- 
long to this group. 

5.3. Comparison between UKIDSS-based and 
2MASS-based disk diagnostics 

Comparing the list of candidate stars with excesses in JHK 
bands selected with £>2mass in GPM07 and the list of stars with 
reliable excesses selected here with Qukidss > we found 

- 139 sources selected only with Q^mass \ 

- 391 sources selected only with Qukidss, among which 307 
are listed in the 2MASS Point Sources Catalog; 

- 184 sources in both lists. 

Among the stars of the first group, 52 have good UKIDSS 
photometry but do not have excesses in Qukidss diagrams. We 
compared the UKIDSS and 2MASS images for each of these 
stars and verified that the excesses found in 2MASS bands for 
all the stars of this group are due not to a disk, but to one or more 
close stars resolved in UKIDSS images but not in 2MASS. As an 
example, Fig. [7] shows one star which belongs to such a group. 
These 52 stars cannot anymore be considered disk-bearing stars 
and therefore cluster members. 

5.4. Contamination of the list of candidate members with disk 

Because not all selected candidate members with disk are X- 
rays sources, a significant contamination by both foreground and 
background sources can be expected. Foreground contaminating 




Fig. 7. UKIDSS (right panel) and 2MASS (left panel) images of 
a disk-less star taken from our multiband catalog corresponding 
to a 2MASS source with K band excess according to 2MASS 
data. 



sources are expected to dominate the sample of 399 stars dis- 
carded from the list of candidate members with disk (see Sect. 
15.21 ). by reason of their optical colors. We retained in our final 
list of candidate stars with disk only ~ 90 stars with infrared ex- 
cesses and optical co lors similar to those of foreground stars. In 
Guarcello et al .1 120101 their nature is discussed in detail, and we 
conclude in that paper that at least 2/3 of these stars are bona-fide 
YSOs associated to Ml 6. 

The identification of candidate background contaminating 
sources is a hard task, because they can be confused with largely 
embedded YSOs associated to M16. However, we do not expect 
a significant amount of background sources among the candi- 
date members with disks detected in the optical bands, because 
the cloud absorbs the background optical radiation efficiently. 
This is evident in the V vs. V - I diagram, where stars in the 
Main Sequence locus are observed up to the distance of M16. In 
order to estimate the amount of visual extinction of background 
sources due to M16, we used the UKIDSS sources without opti- 
cal counterparts, which are dominated by background giant stars, 
because their intrinsic colors can be accurately determined. We 
divided the WFI fields in 400 cells, in which fall between 37 
and 652 sources, accounting for a good statistic. We evaluate the 
median H - K of the sample of sources falling in each cell and 
then the median Eh-k for each cell, taking advantage of the fact 
that the intrinsic H - K color s for giant stars range from 0.07 
to 0.31, with a median of 0.2 dBessell & BretA Il988h . Because 
the sample of UKIDSS sources without optical emission is pop- 
ulated by giant stars a t large distance (up to ~ 10 Kpc, as sug- 
gested by Lucas et ail 120081) . the values of Ay we obtain with 
this procedure is strongly affected by the contribution of the ma- 
terial in the Galactic plane. However, the lowest values we ob- 
tain (8'" < A v < 10'") should to be related to the background 
giant stars which are closest to M16 (and then mostly due to 
the cloud). With this estimate of visual extinction for the back- 
ground stars immediately behind the nebula, we can understand 
whether a background giant star can be observed in our opti- 
cal observation and can contaminate, for some reason, our list 
of candidate members with disk. To this aim, we consider the 
V - K c olor of a GO giant sta r (the brightest type of giant star 
listed in lBessell & Brett|[T988l) . which at a distance of 1800 pc 
and with an extinction ranging from 8 m to 10 m is between 8.35" 1 
and 10.63'". Only 17 candidate stars with disk have the V - K 
color in this range, so we conclude that the contamination by 
background stars of our list of candidate members with disk with 
optical detection is negligible. 

Background giant stars and galaxies observed through the 
Galactic plane can contaminate the sample of candidate stars 
with disk which are not detected in optical, because of the depth 
of the UKIDSS catalog. However, as explained in Sect. 15.21 all 
these stars, except for 21, were classified as "probable cluster 
members" and not used in the study of disk frequency and star- 
formation chronology. In this case, as well, the contamination 
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by background sources is negligible. This is also valid for the 
YSOs selected with the IRAC color-color diagram, as explained 
in GMD09. 

We then can safely assume that our sample of candidate stars 
with disks is not strongly affected by contamination from either 
foreground or background sources. 

6. Final list of cluster members 

We decided to compile the list of stars with excesses in JHK 
including all stars selected with Qukidss with reliable excesses 
(see Sect l5.2l ). plus the stars with only 2MASS excesses with 
saturated UKIDSS photometric quality flags and those with bad 
UKIDSS photometry. Considering also the stars with excesses 
in IRAC bands, candidate stars with disks can be classified as 
follows 

- 627 stars selected with Q indices, among which 

- 269 stars selected with both Qirac an d Qukidss > 

- 74 stars selected only with Qirac', 

- 284 stars selected only with Qukidss ■ 

- 178 stars selected with the IRAC color-color diagram, 
among which 

- 52 stars selected both with the Q indices and the IRAC 
diagram; 

- 126 stars selected only with the IRAC color-color dia- 
gram. 

- 81 stars added from those selected with Q2mass indices, 
which cannot be selected by Qukidss since 

- they are brighter than the saturation limits (5 stars); 

- bad UKIDSS photometric quality flags (76 stars). 

for a total of 834 (627+126+81) candidate cluster members 
with disk (Fig. [8] shows the V vs. V - I and J vs. J - K dia- 
grams with the selected candidate stars with disk). These stars 
with infrared excesses can be divided into three groups. 

- Candidate ClassII objects, from which we observe both the 
emission from the photosphere and from the inner region of 
the disk. Because none among these stars fall in the normal 



stars locus in some of the Q diagram^ the excesses in all the 
NIR bands should come from the same region of the disk, i.e. 
the inner rim at the dusts sublimation radius. 

- Candidate ClassII objects selected only with 2v/y[8.0] (4 
stars), which likely have dissipated their inner disks and their 
Spectral Energy Distributions (SEDs) have a broad silicate 
band in emission, responsible for the excess in [8.0]. 

- Candidate YSOs selected only with the IRAC color-color 
diagram, given the lack of detection in optical or J bands. 
These are likely embedded ClassII and ClassI YSOs. Q in- 
dices are useless to select these stars, while the IRAC color- 
color diagram is an efficient tool to this end. 

The final list of cluster members includes the selected candi- 
date stars with disk (834 stars) and the candidate young stars as- 
sociated with Ml 6 without a disk. These stars are selected from 
the X-ray sources which meet the following requirements: 

- they have at least one stellar counterpart; 

- the positions of the stellar counterparts in the V vs. V — I and 
J vs. / - K diagrams have to be compatible with the cluster 
locus. 

The last requirement discards all the foreground contami- 
nants sources (72), whose position in the color-magnitude dia- 
grams cannot lie in the cluster locus. The first requirement al- 
lows us to discard almost all background contaminating sources, 
whose optical or infrared emission cannot be observed because 
of the absorption by the interstellar material along the galactic 
plane. We estimated the maximum number of expected extra- 
galactic sources observed in our ACIS o bservations taking ad - 
vantage of the logNvs. logS diagrams of Puccetti et al. (2006), 
which is a work based on XMM-Newton observations of the ex- 
tragalactic sources lying in the ELAIS fields. In order to com- 
pare our data with their logN vs. logS diagrams, we convertecQ 
our limit counts rate (6.25 x 10~ 5 cnts/sec) into flux in the 0.5- 
10 keV energy band with a power-law model with energy index 



6 they fall in the loci of stars with excesses or those where excess and 
reddening cannot be discerned. 

7 To this aim we used WEBPIMMS 
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a=0.8 in order to be consistent with lPuccetti et alJ(l2006l) . We set 
the interstellar extinction nH equal to 5.4 x 10 22 cm~ 2 from the 
highest Ay (30'") we can estimate from the sample of UKIDSS 
sources without optical detection, which is dominated by back- 
ground giants stars observed up to large distances (see Sect. 15.4b . 
With this method, we estimate ~ 193.5 extragalactic sources to 
be detected in our ACIS-I observations. But this is clearly an 
upper value to the real number, because we expect further inter- 
vening extinction (i.e. a larger nH) even behind the M16 nebula, 
which reduces the number of detectable extragalactic sources in 
our ACIS-I images, and also for the decrease of limiting sensitiv- 
ity at large off-axis angle in the ACIS-I detector, which has not 
been properly taken into account in the above estimate. Since 
there are 504 X-ray sources without stellar counterparts in our 
catalog, it is possible to assume that about 300 among them are 
very embedded YSOs associated with Ml 6 instead of contami- 
nating background objects. 

The total of candidate disk-less pre-main sequence stars as- 
sociated to the Eagle Nebula is then equal to 1 103 sources. This 
is not a complete census of the X-rays-emitting sources of Ml 6, 
because we lost all quiescent sources with a luminosity Lx below 
our threshold (which is log(Lx) - 29.8 erg/s at 1750 pc) during 
the observation. To estimate the percentage of cluster members 
at various masses that we do not observe, we used the existing 
most complete census of the X-ray population of a 1 Myear clus- 
ter: the Chandra Ultradeep O rion Project (COUP \Getman et al.L 
2005). Fe igelson etaD (120051) calculated the unabsorbed Lx dis- 
tribution in the whole ACIS-I band sampling the Orion sources 
in five different mass bins. With their results and comparing the 
distributions with our detection limit, we evaluated that we ob- 
served 42% of stars with a mass smaller than 0.3 M (lower 
than the limit of our optical and IRAC data), 83% of stars with 
masses between 0.3 and 1 M Q , 92% between 1 and 3 M , 100% 
between 3 and 10 M Q and 75% between the stars more massive 
than 10 M G . We can estimate then that we observed about the 
86% of clusters members more massive that 0.3 M Q , which is 
the limit of our optical and IRAC data (and which is the limit 
of the members selected to study the spatial variation of disks 
frequency and the chronology of star formation in the following 
sections). 

7. Spatial variation of disk frequency 

Disk frequency is expected to vary across the Eagle Nebula given 
the presence of recent star-formation sites and the not uniform 
spatial distribution of massive stars, which are mostly concen- 
trated in the central cavity corresponding to NGC 6611. The 
spatial variation of disk frequency in the central region (the C- 
field) has been studied in GPM09 and GMD09, where it was 
found that the disk frequency decreases close to massive stars. 
These studies pertained only to the C- field, which contains al- 
most all the massive stars of NGC 661 1. In this field, we selected 
here 896 candidate PMS members without disk and 510 candi- 
date cluster members with disk. The new list of cluster members 
compiled in this work with the UKIDSS data confirms the de- 
crease of the disk frequency close to massive stars, as shown in 
the histogram (left panel of Fig. |9j of the disk frequency vs. the 
flux emitted by OB stars and incident on cluster members. The 
incident fluxes (see GPM07 for detail on flux evaluation) were 
calculated by summing the emission from all massive stars in 
the cluster and incident on each cluster members brighter than 
J = IT" (i.e. with masses larger than 0.1 - 0.2 M , in order 
to explore the same mass ranges of disk-less and disk-bearing 
candidate members). The dotted lines in the left panel of Fig. 



[9] show the disk frequencies in each UV flux bin obtained in 
GMD09. The disk frequencies found in this work are higher by 
~10% than those found in GMD09. This reflects the larger num- 
ber of candidate disk-bearing stars we selected with UKIDSS. 
The average disk frequency in the whole C- field is now equal to 
36% +1% (33% ± 2% taking into account the completeness of 
our X-ray survey for stars more massive than 0.3 M , see Sect. 

The visual extinction map, discussed in AppendixlBl allowed 
us to obtain the dereddened Vo vs (V - 7) and / vs (J - K) di- 
agrams. Thanks to these diagrams, we e stimated th e mass es of 
the cluster members from the tracks of Siess et al. (2000) and 
we studied the disk frequency at different mass ranges. Table [3] 
shows the average disk frequency in the C-field (which con- 
tains NGC 6611) for two different mass ranges. There is no 
large dependence of the disk frequency on the stellar mass, 
even if a slightly (1.5 x (Tf req ) higher frequency is found for 
masses < 1M . This difference becomes even less significant 
taking into account the percentage of X-ray sources below our 
Lx threshold (Sect. [6). The right panel in Fig.[9]shows the disk 
frequency for two different ranges of central star masses (larger 
and smaller than 1M , respectively) and for different values of 
incident fluxes emitted by massive stars, taking into account only 
the candidate members that fall inside the C-field. The differ- 
ence between the disk frequencies for high- and low-mass stars 
does not strongly depend on the intensity of the incident flux. 
Note that the disk frequencies in Table [3]differ slightly from the 
average frequency in NGC 661 1 (36% + 1%), because it was im- 
possible to evaluate the masses of all candidate cluster members 
that fall in the C-field. 

Table 3. Disk frequency in the C-field for different mass ranges 



Mass range 


Disk-less 


Disk-bearing 


Disk frequency 


> 1M 


119 


56 


32% ± 4% 


0.3-1 M Q 


394 


300 


43% ± 3% 



In order to compare with a good statistic the spatial variation 
of disk frequency in the outer regions of M16 with the morphol- 
ogy of the nebula, the average age of cluster members, and the 
positions of OB stars, we defined 20 regions across the portion of 
the Eagle Nebula observed with ACIS-I. Figure \W\ shows these 
regions and their names, the disk frequencies (on a gray scale) 
and their content of OB stars. Each region is labeled with the 
name of the ACIS-I field in which it is contained (the first let- 
ters, one or two), and the orientation with respect to the center 
of this ACIS-I field (for instance, the region NE_S is the region 
south of the NE- field). Figure[10]shows that the disk frequency 
is higher in specific outer regions of the nebula with respect to 
NGC 6611 (mostly contained in the C-C field). This could be 
due either to the presence of a younger population of members 
or to a lack of OB stars in the outer regions. The disk frequency 
in all regions varies from 52% + 13% (C.SE region) to 17% + 8% 
(E_SE region). 

Finally, as explained in Sect. 15.21 we used Ry = 3.1 in the 
definitions of Q indices. However, it is well known that the red- 
dening law in the direction of NGC 661 1 is anomalous, even if 
GPM07 revised the estimate of Ry, finding a value (3.27) only 
slightly different from the canonical one. Performing the selec- 
tion of disks with Ry = 3.7, the number of selected candidate 
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Fig. 9. Left panel: histogram of the disk frequency in NGC 6611 (the C- field) vs. the flux emitted by massive stars and incident 
on cluster members, compared with that found in GMD09 (dotted histogram). Right panel: histograms analogous to that in the left 
panel, but with the cluster members in two different ranges of stellar masses: > IM Q (solid line) and O.2M < M stars < 1M (dotted 
line). 



members with disks decreases to 8%, and the average disk frac- 
tion to 3%. The spatial variation of the disk frequency and its 
relation with stellar masses remains essentially unchanged. 

8. Chronology of star formation 

In order to shed some light on the star-formation history in the 
Eagle Nebula and to understand if the disk frequency correlates 
with the age of the stellar population, we e valuated the age o f 
each disk-less member from the isochrones of Siess et al. (2000), 
with the Vq vs (V - 7) and Jo vs (J - K) diagrams. We ex- 
cluded disk-bearing members from this computation, because 
their colors were not purely photospheric. A total of 333 disk- 
less members can be plotted in both diagrams, obtaining two 
different age estimates for these stars. On average, optical and 
infrared age estimates for each star agree, with a median differ- 
ence of ~ 0.02 Myears but with a large spread. The age of these 
stars was then set equal to the average value of the two estimates. 

FigureQT]shows the dereddened Vo vs. (V - T)q diagrams for 
the stars falling only in the C- field (upper left panel); in the sky 



region common to C- field and E- field (upper right panel); only 
in the E-field (lower left panel); and, finally, for the sources 
falling in the NE-field. The sequence upper-left, upper-right, 
lower-left panels go through the nebula from NGC 66 1 1 east- 
ward. In all these diagrams, crosses mark the X-ray sources. 
The disk-less members, traced by X-ray active stars, in the upper 
right panel (region in common between the two ACIS-I fields) 
are apparently older than those in the C- field, and younger than 
those in the E-field. The median age for these three groups of 
stars, indeed, increases toward the east: 1 .0 Myear for the stars in 
the C- field; 1.4 for the stars in the common region; 2.0 Myears 
for the stars in the E-field. This trend contradicts the hypoth- 
esis of a large-scale triggering of star formation in Ml 6 by the 
massive stars in its center, because then the oldest stars should 
be found in the C- field. Stars in the NE-field are more sparse 
in the diagram, with a fraction apparently younger than 1 Myear. 
The Jo vs. {J-K)o diagrams of the same regions, not shown here, 
share the same properties as the diagrams in Fig.QT] Following 
the extinction map shown in Appendix [B] Ay changes slightly 
going from the center of the field eastward, with values between 
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Fig. 10. Gray-scale map of the disk frequency in M16. Dots and 
stars mark the positions of OB (earlier than B5) stars. The dotted 
boxes are the ACIS-I fields. The circles mark the regions we 
defined. 



Stars only in the C-field Stars both in C- and K-fields 




Fig. 11. Intrinsic Vq vs. (V - 7)o diagrams for the optical sources 
(points) and the optical+X-ray sources (crosses) that fall in four 
different sky regions: only in the C-field (upper left panel), in 
the common region of C-field and E- field (upper right panel), 
only in the E- field (lower left panel) and fall ing in the NE- fiel d 
(lower right panel). The isochrones are from Siess et al. (2000). 



2 m -4"'. In particular, it is on average smaller in the east direction 
with respect to the center. This allows us to exclude the hypothe- 
sis that this result can be influenced somehow by contamination 
of the member list (which, however, is not crucial, as explained 
in Sect. [6). Indeed, if most of the sources in these fields were 
background sources, their extinction (and then their age) would 
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Fig. 12. Map of the age of disk-less cluster members. The 
positions of Classl/embedded ClassII YSOs are also shown 
(ClassNNIDs are sources whose IRAC colors allow both clas- 
sifications). For each region we give the median age of disk-less 
sources and the number of stars considered for age determina- 
tion. North is up, East is left. 

be underestimated, which mostly affect the sources in the East, 
which are older that those in the central field in Fig.fTTl 

In order to investigate the triggering of star formation in Ml 6 
by the massive stars associated with NGC 661 1 on a large spatial 
scale, we constructed a map of the age of disk-less cluster mem- 
bers in different regions of M16 dividing the WFI FoV in 25 cells 
of 1 .7 x 1 .4 parsecs. The representative age of each cell was cho- 
sen as the median age of all disk-less members falling in that 
cell. The cell in which most of NGC 6611 falls is 0.8 Myears 
old, younger than all surrounding cells in every direction but 
northward. Among the cells with a significant number of stars, 
the oldest cells (2-3 Myears old) are in the East and South; the 
youngest (0.1 Myears) is northward from NGC 6611. This grid 
is not shown here. However, despite an in-outward sequential 
star formation in M16, which should to be observed in the case 
of a large-scale triggering of star formation by the massive stars 
associated to NGC 6611, the grid we obtained suggests an age 
trend from southeast to northwest. This age gradient is shown 
in Fig. [12] This map contradicts the hypothesis of a large-scale 
triggering of star formation by the massive stars of NGC 6611. 
The oldest disk-less cluster members are in the southeast re- 
gion (2.6 Myear), while the youngest are in the northwest region 
(0.3 Myears). Stars in the region of M16 which corresponds to 
NGC 6611 are 1.3-0.8 Myears old. This trend is confirmed not 
only by the median values, but also by the "global offset" of the 
whole age distributions toward younger ages going from south- 
west to northeast. Figure [12] also shows the spatial distribution 
of Classl/embedded ClassII objects. These stars are sparse, with 
a group clustered in the northeast, in the center of the NE-field. 
This group corresponds to the position of the embedded NE clus- 
ter (see Sect. [TJ and the chronology of star-formation we found 
suggests that this is a young star formation site in M16. 

9. Discussion 

9.1. Disk frequency and externally induced photoevaporation 

In Sect. [7] we found an average disk frequency in the C-field 
equal to 36% + 1%. In the same sky region, GMD09 evaluated 
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an average disk frequency of 24% + 2%, lower than our value. 
This discrepancy is not surprising, since the new IR photometry 
is deeper (about 1 .5 m in J) and the new list of members is larger 
than that found in GMD09, in agreement with the smaller pho- 
tometric errors in UKIDSS and the higher sensitivity and spatial 
resolution in UKIDSS/GPS with respect to 2MASS/PSC (allow- 
ing a more complete selection of disk-bearing stars). Besides, 
as explained in the next section, the disk frequency is expected 
to increase for stars with decreasing masses. In GMD09, the 
obtained disk frequency in the central re gion of NGC 6611 
was compared with the value estimated by lOliveira et al. 1 120051 
(~ 58%), with the conclusion that the discrepancy is due to the 
deeper infrared observation used in the latter work, and not to an 
inefficiency of the Q method to select stars with a disk. The new 
average disk frequency estimated in this work, which is higher 
than that found in GMD09, supports this hypothesis. 

In order to search for a possible relation between disk fre- 
quency and stellar mass, we evaluated the mass of the can- 
didate cluster members from the de-reddened Vo vs. (V - I)o 
and Jo vs. (J - K\) diagrams (see Sect. |7). We limited this 
analysis to the C- field because NGC 6611 is completely con- 
tained in this field. Table [3] shows the average disk frequen- 
cies in two bins of stellar mass. A higher frequency has been 
found for stars with mass between 0.2 and 1 M Q , in agree- 
ment with the results ob t ained in other star- forming regions , 
as 77 Cha dMegeath et all 120051). Upper S e p dCarpenter et all 
2006) , Trumpler 37 (ISicilia-Aguilar etafl. 120061) and IC348 
dLada et alll2006l) . However, these studies evaluated for the pop- 
ulation of stars more massive than IM Q disk frequencies from 
few percents to ~ 10%, which is significantly lower than our 
38%. This difference is in part explained by the different meth- 
ods used to select disk-less and disk-bearing members. Most of 
this discrepancy, however, is due to the different age of the clus- 
ters. Indeed, all these clusters have an age ranging between 3-9 
Myears, older than the members of NGC 6611 (whose median 
age is lower than 1 Myear). The disk frequency found in our 
work consequently refers to a younger pre-main sequence popu- 
lation, with a larger number of intermediate- to high-mass stars 
that have not yet dissipated their disks. 

With the more complete disk selection of this work, the trend 
of disk frequency with the incident flux found in GPM07 and 
GMD09 changed slightly (as shown in the left panel of Fig. |9). 
In that work the disk frequency decreases monotonically from 
larger to smaller distances from OB stars. With the new and 
more reliable list of cluster members, the disk frequency is al- 
most constant in the field, dropping only close to the massive 
stars, where disks are irradiated by intense UV radiation. Our 
disk selection is sensitive to stars with an inner disk and those 
with a population of small grains in their disk. The decrease of 
the percentage of disks selected by NIR excesses means a dis- 
sipation of the whole disks (down to the inner region) and/or of 
the small (~ some fori) grain population. iGuarcello et al.l (12009) 
discussed several processes which lead to a decrease of the small 
grains population in photoevaporating disks. 

When externally induced photoevaporation occurs, disks are 
destroyed from outside (where the disk surface density is lower) 
inward, down to the gravitational radius (R e ) where the es cape 
velocity is equal to the sound speed (see Hollenbach 1994). For 
instance, in a disk around a 0.2 M Q star heated up to 1000 K 
by incident Far-Ultraviolet (FUV) radiation, R g is at 20 AU. 
R g decreases significantly when the disk is irradiated by the 
Extreme Ultraviolet (EUV) radiation and it is heated up to 
~ 10000 K. The dissipation of the inner disk is also more rapid 
in EUV regime, occurring in timescales smaller than 1 Myear 
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Fig. 13. Disk frequencies (in percenta ge) vs. age for the stellar 
population of the clusters studied in lHaisch et alj|200ll (black 
filled dots) and the regions of the Eagle Nebula defined in Fig. 
[TUl (empty circles). The line marks the age vs. disk frequency 
relation interpolated from the Haisch data. The size of the empty 
circles is related to the content of OB stars in the regions, as 
explained in the figure. 



( Storz er & Hollenbachl 1 19991) . This is crucial in our study, be- 
cause our disk diagnostic is sensitive to the emission from the 
inner region of the disk and the average age of stars in NGC66 1 1 
is ~ 1 Myear. Extreme ultraviolet radiation dominates the heat- 
ing of the disks close to OB stars, up to an average distance of 
about 1 parsec. The left panel of Fig. [9] shows that the disk fre- 
quency drops close to OB stars, while it is almost constant at 
larger distance. This is consistent with the observation that in 
NGC 6611, with an age of ~1 Myear, only disks irradiated by 
EUV radiation have been significantly dissipated down to the 
inner region, while FUV radiation in the outer part of the cluster 
ha s not affecte d the d isk evolution. Our results agree with those 
of iBalog et all (120071) . who found that in NGC 2244 (a young 
cluster hosting a O star) the disk frequency drops at a distance 
smaller than 0.5 parsec from this O star, and with the HST ob- 
servations in the Trapezium in ONC, where photoevaporating 
disks have been observed at distance from the O stars smaller 
than 0.4 parsec. The bin with the highest UV flux in Fig.[9]in- 
deed, corresponds to an average distance from massive stars of 
~1 parsec. These results go in the direction that only EUV ra- 
diation can affect the planetary formation in circumstellar disks, 
because photoevaporation in the FUV regime could be important 
only on timescales larger than the lifetime of O stars or than the 
cluster relaxation time. 

We also searched for a relation between the observed disk 
frequency and the median age of cluster members in the regions 
of the Eagle Nebula defined in Fig. [TOj To this end, we com- 
pared the observed disk frequencies with th ose expected on th e 
basis of the age of the clusters analyzed in lHaisch et a l. (2001). 
Figure [T3lshows the d isk frequencies of both the clusters studied 
in Haisch et al.ll200ll (filled circles) and the region of the Eagle 
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Nebula defined in Fig. [10] (empty circles). The size of the empty 
circle is related with the amount of massive stars in the region, 
as shown in the figure. The disk frequency of these regions can 
be easily compared with the interpolated age vs. disk frequency 
relation, which is marked by the line. Because three of the four 
regions with the difference between the observed frequency and 
the corresponding expected frequency (from the interpolated re- 
lation) larger than 60% are among those with the largest popu- 
lation of OB stars, a connection between the number of massive 
stars in the sub-fields and the difference between "expected" and 
observed disk frequencies can be assumed. 

Note that the regions with the highest differences are also the 
youngest. It is possible that in these regions most of the associ- 
ated YSOs with disk are very embedded in the cluster medium, 
which makes their detection hard and partially accounts for the 
observed differences. The only field in which this hypothesis 
does not work is the C-C field, which corresponds to the cavity 
cleared by NGC 6611 stars (and among the region with lower 
extinction, see Appendix IB1. In C-C field the difference is more 
likely due to the high number of massive stars. In conclusion, 
Fig. [13] shows a larger difference between the observed and the 
"expected" disk frequency for those regions that are younger av- 
erage and for those with more massive stars. 



9.2. Efficiency of disk photoevaporation for different masses 
of the central star 



The right panel of Fig.[9]shows how the disk frequency for low- 
and intermediate- to high-mass stars varies with the incident flux 
emitted by OB stars. This analysis is performed to study the rel- 
ative role of two different mechanisms for disk dissipation which 
can take place in young massive clusters. For a unperturbed evo- 
lution, the timescale for disk dissipation strongly depends on the 
mass of the central star. Even if there is not yet an universal 
scaling law which gives the time of dissipation of the inner disk 
versus the mass of the central star, high-mass stars emit more en- 
ergetic flux and more intense stellar winds than low-mass stars, 
therefore they erode their inner disks faster. Indeed, the disk fre- 
quencies for high-mass stars decline faster than those evaluated 
on the entire cluste r populations (which declines in ~ lOMyears, 
Haisch et al. 2001). On the other hand, for disks photoevapo- 
rated by close massive stars, specific simulations predict a faster 
dissipation of disks in low-mass stars w ith respect to high-mass 
stars. For instance, Adams et al. (2004) have shown that the in- 
tensity of the flux of photoevaporating gas from disks with a cen- 
tral star of O.25M is less intense (by one order of magnitude) 
than those with a star of 1M . This is due to the weaker gravita- 
tional field produced by less massive stars, which also results in 
a larger gravitational radius. 

In the right panel of Fig. [9] the difference between the fre- 
quencies for high- and low-mass stars is almost independent 
from the intensity of the incident flux. This suggests that in 
NGC 661 1 the disk dissipation timescales due to photoevapora- 
tion induced by OB stars do not strongly depend on the mass of 
the central stars. This can be explained by the age of the cluster 
(~1 M year). Indeed, following the simulations of I Adams et al.l 
(2004)' and using the typical mass-loss rates which they evalu- 
ated, the total mass lost in 1 Myear from two disks around a 1 
M Q and 0.5 M star, respectively, differs by less than ~ O.O1M . 
This low value cannot produce observable differences in NIR 
excesses. 
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Fig. 14. DSS-I image of about 4° x 3° of the sky region between 
Ml 6 and M17. The axis of the galactic coordinates, the limit of 
the molecular shell (black curve), and the limits of the regions in 
Fig. [12] (white lines) are shown. 



9.3. Large-scale star-formation triggering induced by OB 
stars 

In Sect. [8] we tested the hypothesis of a possible large-scale 
triggering of star-formation in the whole M16 by the massive 
stars in NGC 661 1, searching evidence that younger objects are 
concentrated in the outer regions of the Eagle Nebula. The age 
trend shown in Fig.[T2]negates this hypothesis: only in the north- 
ward direction the disk-less stars associated with the nebula are 
younger than those in the central region, while in the other di- 
rection older objects have been found. Instead, there is a gra- 
dient trend in the SE-NW direction, with a difference of 2.3 
Myears between the two extremes (see Fig.[T2l. Hence, the pic- 
ture which arises from this is that the star-formation events in the 
southeast region of M16 are older (of 1-2 Myears) than those in 
NGC 661 1 (central region). Only the star formation in the North 
and West (i.e. the region around the massive star in formation 
and the SFO30 cloud) may have been triggered by the radiation 
from the massive stars of NGC 6611 in recent time (less than 
1 Myear), as expected by a small-scale triggering phenomenon. 
The star formation in the embedded cluster in the North-East 
occurred in the last 0.5-0.3 Myears. 

If this sequence of star formation is correct, something has 
given the initial impulse for the star formation in the South-East 
about 3 million years ago. We tentatively identified this event 
in the incidence on M16 of a giant molecular shell discovered 
bv lMoriguchi et akl (120021) . Following these authors, this shell is 
about 170 x 250 parsecs wide, and was created about 6 Myears 
ago likely by supernova explosions. This shell reached the Eagle 
Nebula in the southern direction approximately 2-3 Myears ago. 

Figure [14] shows a DSS-I 4° x 3° image of the region of 
the galactic plane with the Eagle Nebula (Ml 6) and the Omega 
Nebula (Ml 7). The limits of the molecular super-shell discov- 
ered bv lMoriguchi et alJ 2002 (the black curve) and the limits of 
the diagonal regions in Fig.[T2](the white lines) are also shown, 
with the oldest region of Fig. Q~2] closer to the shell. Both the di- 
rection and the time at which the molecular shell encountered the 
Eagle Nebula are consistent with the sequence of star-formation 
suggested by this work. Moreover, both the star-formation events 
in the southeast region of M16 and the oldest events in the 
Omega Nebul a (M17) can be dat ed to approximately 3 million 
years ago (see lPovich et all 12009). which supports a connection 
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between the stars-formation events in the two star-forming re- 
gions. 



10. Conclusions 

We studied the chronology of star formation in the Eagle Nebula 
(Ml 6) and the spatial variation of disk frequency to test the ef- 
fects of OB stars on the parental cloud and on the nearby pop- 
ulation of young stars. We analyzed optical, infrared, and X-ray 
data of this region. X-ray data have been obtained from the ACIS 
archive (one observation centered on NGC 6611 in the center 
of Ml 6) and from two new observations centered on the pillar 
structure called "Column V" and on a recently discovered em- 
bedded cluster (in the East and northeast regions of Ml 6, re- 
spectively). We detected 1755 X-ray sources inside these three 
regions. Infrared data in JHK bands of the whole region have 
been obtained from the UKIDSS/GPS catalog. Optical BVI and 
Spitzer/IRAC data of Ml 6, used in this work, have been pre- 
viously analyzed in other publications. With all these data, we 
compiled a multiband catalog of this region. 

We selected the candidate stars associated with M16 using 
two criteria: infrared excesses for disk-bearing sources and X- 
ray emission for disk-less stars. Infrared excesses have been de- 
tected with two different diagnostics: the IRAC color-color di- 
agram and optical-infrared color indices defined in order to be 
independent from the extinction. We selected a total of 834 can- 
didate disk-bearing stars. The list of young stars associated with 
M16 has been completed with 1103 candidate disk-less mem- 
bers, for a total of 1937 sources. 

We have found an average disk frequency of 36% + 1 % in the 
central region of the nebula, which corresponds to NGC 66 1 1 . 
We also found a weak dependence of the disk frequency on stel- 
lar mass: the disk frequency is slightly higher for stars less mas- 
sive than 1M Q with respect to more massive stars. For stars more 
massive than 1 M Q , the frequency is equal to 32% + 4%, higher 
than the values observed in older star-forming regions, likely be- 
cause in 1 Myear (which is the median age of NGC 661 1) a sig- 
nificant number of high-mass stars have not yet dissipated their 
disks. 

Disk frequency drops at small distance (d< 1 pc) from the 
massive stars of NGC 6611, where disks are irradiated by in- 
tense EUV fluxes. The variation of the disk frequency with the 
intensity of incident UV flux does not depend on the mass of 
the central stars. These results indicate that in ~1 Myear photoe- 
vaporation induced externally by EUV radiation has completely 
dissipated only disks close (d < 1 parsec) to OB stars, with- 
out observable effects at larger distance (where photoevaporation 
is driven by FUV photons) and without an evident dependence 
of the efficiency of photoevaporation on the mass of the central 
stars. This study supports the hypothesis that only EUV radia- 
tion from OB stars may have an effect on the planet formation 
inside circumstellar disks. 

We studied the history of star formation in the nebula, finding 
that YSOs in NGC 66 1 1 are on average younger than those in the 
regions of M16 in the East and South, and older than those in the 
North. We derived a chronology of star formation in the nebula 
from southeast (about 2.6 Myears ago) to northwest (less than 1 
Myear ago). This chronology of star formation does not support 
the hypothesis that OB stars in NGC 6611 have triggered the 
formation of new stars in the whole nebula up to a distance of 
~10 parsec. We explained this chronology with a star formation 
process induced externally, and we tentatively identified a giant 
molecular shell, which reached M16 about 2-3 million years ago, 



starting the star formation, and which could connect the Eagle 
Nebula with the Omega Nebula. 
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Table A.l. Symbols used to mark the stars in Figs. lA.lllA.21 and lA.3l 



Type of star 



Symbol 



X-ray sources with a stellar counterpart X 

X-ray sources without any counterpart pluses 

Stars with disk added from 2MASS member list circles 

Stars with excesses only in UKIDSS bands squares 

Stars with excesses in IRAC bands diamonds 




Fig. A.l. VIMOS@ESO image of the ColumnV in R band. Symbols are described in Table IA. II 



Appendix A: YSO population in selected regions of M16 

In this appendix we briefly describe the candidate members of Ml 6 associated with the ColumnV, with the embedded cluster in 
the northeast and with the cavity at northwest. In the images of these regions, Figs. IA.lllA.2l and lA.3l the candidate members are 
marked with the symbols described in Table lA.fl 

Figure [AH shows a 10 sec image in R band of the ColumnV, obtained with VIMOS@ESO on 2009-03-24 (obs. ID: 083.C- 
0837, PI: Guarcello) as part of the VIMOS pre - imaging of the Eagle Nebula. The sky region in this image is 7.8' x 3.6' wide 
(corresponding to 3.7 x 1.7 parsec), and it is 7' (3.3 parsec) eastward from NGC 6611. In this region are 97 X-ray sources (among 
which 23 have also infrared excesses, 22 do not have any stellar counterpart, 52 are candidate disk-less YSOs), 52 candidate disk- 
bearing stars (31 with only UKIDSS excesses; 18 with IRAC excesses; 3 added from the 2MASS list of members), for a total of 104 
candidate members and an average disk frequency equal to 50% ± 8%. In the various color-color and color-magnitude diagrams, the 
positions of all sources in Fig. lA.ll indicate that they are on average affected by the typical visual extinction of the optical sources 
associated with the nebula (Ay ~ 2.5'"-3'"), and that they are on average 1-1.5 Myears old. Following the star-formation chronology 
across the Eagle Nebula shown in Fig. [12] it can be assumed that the ColumnV and the Elephant Trunks have been formed almost 
simultaneously in the last 1-1.5 million years. Indeed, the few sources younger than 1 Myears are those in front of the cap of the 
pillars. Th e massive stars closest to the ColumnV are W541 (a =18:19:19.123, 6 =-13:43:52.32), a binary dMartavan et al.Ll2008[) 
KL5 star dHillenbrand et al.L [19931) 1.3' northward, and W472 (a =18:19:04.709, 6 =-13:44:44.52), a B3 star dHillenbrand et all 
1993), in correspondence of the cap of the pillar. The source W541 may be responsible for the erosion of the center of the pillars, 
while W472 may be the ionizing source which heats the cap of the pillar and may have induced the m ore recent star-formation 
events. In Fig. IA.ll the letter "A" marks the position of four water masers identified bv lHealv et al . (2004). One embedded candidate 
ClassI YSO, with very red [3.6] - [4.5] color, can be associated with these water masers (see also Ind ebetouw et all 120071) . We also 
identified an X-ray source without stellar counterpart, 1.6" - 2" away from the water masers. This source could also correspond to 
a very embedded YSO associated with the water masers, also because of its hard spectrum (its median photon energy is equal to 3.5 
keV). The region in front of the pillar is richly populated: 7 disk-less X-ray sources and 14 candidate sta rs with disk (among which 
8 X-ray sources). The point "B" in Fig. lA.ll marks an ionized knot candidate to be a Herbig-Haro object (Mea burn & WalsnU l986). 
A soft (median energy equal to 1.9 keV) X-ray source and a faint stellar counterpart is 1" northward from the knot. 

Figure IA~2l shows the IRAC [8.0] image of the sky region corresponding to the embedded cluster in the North-East (7.3' x 5.3' 
wide, corresponding to 3.4 x 2.5 parsec). This region includes 59 X-ray sources (among which 24 have no stellar counterpart, 9 
have infrared excesses and 26 are disk-less members), and 25 candidate stars with disk (21 with excesses in IRAC bands, 4 only 
in UKIDSS bands). Both the extinction map, the color-magnitude diagrams and the age trend (Figg. IB~2l and [T2l show that this 
region is populated by very young (age < 1 Myear) and embedded sources (4™ < Ay < 2Q m ). Disk-bearing stars in this region are 
heavily extincted, so they lack optical counterpart (except one star with faint / emission) and they are very red in UKIDSS bands 
(J - K > 3 m ). Also in the IRAC color-color plane, objects in this region can be classified as very embedded ClassII/ClassI sources 
(14 stars), with also 4 ClassII sources with moderate reddening and 1 ClassI YSOs. A very hard X-ray source is associated with this 
ClassI YSO, with 31 net counts and a median energy of 5.5 keV. According to the high extinction, in this region the X-ray sources 
with infrared counterpart have very red colors and, on average, hard spectra, with typical median energies of 3 keV. 
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Fig. A.2. IRAC [8.0] image of the embedded cluster in the North-East. Symbols are described in Table lA.il 




Fig. A.3. IRAC [8.0] image of the region in the northwest of the Eagle Nebula. Symbols are described in Table I A. fl Boxes are the 
ACIS-I FoVs, while the dashed line marks the sky region inside the cavity, close to W584, where no disk-bearing stars have been 
selected. 

Figure lA3l shows the sky region in the northwest of NGC 661 1 with the side of approximately 10' (4.8 parsecs) which has not 
been observed in X-ray (the solid lines in Fig. lA.3l indicate part of the ACIS-I FoVs). In this region, M16 has a cavity-like structure, 
with the 09V W584 star near the center of this cavity, at ~ 0.7 parsec from the well defined NW border of the cavity. This kind of 
structure could have been created by the intense wind emitted by this massive star. As shown in Fig. [A3] no star with disk falls close 
to W584, while 20 candidate disk-bearing members are selected in every direction along the border of the cavity. The formation 
of these stars could have been induced by the compression of the nebula by W584. Since it is more likely that the nebula has been 
swept away continuously from the position of W584 outward, formation of new stars should be induced also inside the cavity, and 
not only in its present borders. Then, it is possible to suppose that W584 has not only created the cavity, but also induced a fast 
dissipation of the circumstellar disks in these young stars. An X-ray follow-up could be useful to select candidate disk-less YSOs 
inside the cavity, to clarify this point. 
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Fig. B.l. Spatial distribution of the UKIDSS sources (gray points) and of the X-ray sources (black dots). Dashed boxes mark the 
ACIS-I fields. North is up, East is on the left 



Appendix B: Extinction map of the Eagle Nebula 

In Sects. [7] and [8] we used the individual extinction of cluster members in order to derive both stellar masses and ages. It was not 
possible to adopt an average extinction for all stars associated with M16 because the density of the cloud varies significantly across 
the field, as it is evident from the spatial distribution of UKIDSS sources (marked with gray points) shown in Fig. IB. II Then, stellar 
extinction values have been evaluated by means of a map of visual extinction Ay toward the Eagle Nebula. Figure IbTTI also shows 
the spatial variation of X-ray sources (black dots). Because they are dominated by cluster members, X-ray sources are clustered on 
NGC 661 1 (in the center of the C- field) with a smooth decrease of source density eastward. The region in the NE- field, where falls 
a group of X-ray sources in a region of high nebulosity, corresponds to the site of the northeast embedded cluster (see AppendixlAli. 

Usually, using photometric data to derive the extinction map of a given region, it is necessary to identify a sample of stars whose 
intrinsic colors are accurately known. In this way, the interstellar reddening can be evaluated as the difference between the observed 
and t he intrinsic colors of these stars. Usually giant stars, whose average infrared colors have been evaluated by several authors 
(i.e. iBessell & Bret! 1 19881) . are used to do this. In our case, giant stars dominate the sample of UKIDSS sources without optical 
counterpart, but their extinctions are dominated by the unk nown distribution o f background interstellar medium in our Galaxy from 
M16 up to very large distances (~ lOKpc, as suggested by iLucas et al.L [20081 and by the J - H vs. H - K diagram for stars in the 
WFI FoV), and not by the material associated with M16. 

It is possible to obtain an accurate estimate of the visual extinction toward the Eagle Nebula by considering the disk-less X-ray 
sources with NIR emission, because they are mostly PMS stars associated with the nebula and their emission is not affected by 
the background inte rstellar medium. Typically, T-Tauri stars with an age of 1 Myear and H - K > 0.2 m are of K-M s pectral type 
(lMeveretal.Lll997l) . Their intrinsic ( J - H) ranges from 0.5 m to 0.7 m , with a median value of 0.66 m (ISiess et al.ll2000l) . Assuming 
this value, it is possible to estimate the individual Ay from the relation Ej-h — (J - H)observed - (J - H)o and Ej-h = 0.107 x Ay. In 
order to obtain reliable results, it is mandatory to exclude candidate disk-bearing stars, whose infrared colors are not photospheric. 
Our selection of disk-bearing stars is discussed in detail in Sect. [5] Starting from the individual Ay, we derived the extinction map 
by dividing the WFI FoV in 100 squared cells and we assigned to each cell the median value of the extinction values of all selected 
stars falling in that cell. 

Figure lB~2l shows the M16 extinction map obtained with the technique described above. The extinction map shows an increase 
of visual extinction from the central cavity, corresponding to NGC 661 1, northward. Ay peaks in correspondence of some structures 
of the nebula describ ed in Sect. [2] as the cluster in the North-East. The increase of extinction northward has been already found by 
iBelikov et al.l d 1999b . 

With the standard hydrogen column (Nh) vs. visual extinction relation dPredehl & Schmittlll995l) . the difference of the average 
extinction between the central and the outer regions of the nebula implies an average hydrogen column density which ranges from 
4.7 x 10 21 atoms/cm 2 in the center to 3.2 x 10 22 atoms/cm 2 in the denser regions. 

In order to verify that the drop of the disk frequency close to the massive stars (see Fig. |9]l is not due to a larger numbers 
of background contaminants selected as disk-bearing cluster members in the outer region of NGC 6611, resulting in higher disk 
frequency values, we compare the visual extinction inside the C - field, where NGC 661 1 lies, with the spatial variation of the UV 
field using Fig. IB. 31 The left panel shows the Ay map, which is a magnification of the central part of Fig. IB. 21 The gray-scale is 
also the same. The right map has been created evaluating the UV field incident on the position corresponding to the center of each 
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Fig. B.2. Ay map in the WFI FoV obtained from candidate disk-less members of M16. The dotted boxes are the ACIS-I FoVs 
analyzed in this work. Ay cannot be computed in the white boxes because of the lack of disk-less members with adequate quality 
photometry in J and H bands. North is up, East is on the left. 




Fig. B.3. Left panel: map of the visual extinction of the central region of M16, corresponding to NGC 661 1. The gray scale is the 
same of that used in Fig. IB. 21 Right panel: map of the intensity of the UV field across NGC 661 1. The gray scale relates the median 
intensity of the UV field to the bins in the left panel of Fig. [9] In both panels the inclined box delimits the ACIS-I FoV. North is up, 
East is on the left. 



bin, taking into account the flux emitted by all massive stars of NGC 661 1. The gray scale relates the intensity value to the bins in 
the histogram in Fig. [9] Comparing both images, no correlation between the intensity of the UV field and the visual extinction, i.e.: 
with the density of the cloud, arises. This leads to the conclusion that our result about the spatial variation of disk frequency is not 
affected by background contamination. 
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Appendix C: Cross-correlation procedure 

Here we describe the procedure we adopted to cross-correlate all catalogs in the field of Ml 6 used in this work. The match is 
performed by identifying the sources belonging to different catalogs with a separation smaller than a given matching radius. The 
used matching radii have been chosen as the highest values at which the spurious matches are negligible with respect to the real 
cross-identific ations. The expecte d number of spurious matches at a given matching radius have been evaluated with the method 
described in Damiani et al 

In order to build the joined catalog of X-ray sources, we noted that the overlap region among E- field and NE- field is very small 
(Fig. IB. 1 1 ), with only few sources found in both fields. Given the extent of the overlapping region between the C- field and E- field, 
we found the common sources by matching the corresponding sources list, adopting a matching radius equal to 1" for sources with 
off-axis smaller than 7' and 2" otherwise (considering the highest value of the source off-axis angle). A total of 76 sources was 
found in both fields, with no multiple matches. To match the optical and the infrared catalogs we used a matching radius of 0.3". 

Because in the ACIS-I field the PSF degrades at large off-axis, to match optical+infrared and X-ray catalogs it was mandatory to 
define a matching radius which depends on the sources off-axis (0) in the X-ray images. Because the uncertainty on source positions 
grows with the off-axis, as a consequence of the spreading of the PSF, we defined the matching radius (in arcseconds) as 

r matc h=Axf(6), (C.l) 

where f(G) is an interpolated relation for the source's off-axis and the position error derived by PWDetect. Equation |CT| holds 
for off-axis> 3', because for smaller off-axis the position errors do not vary with 6. To find the value of A, we should compare the 
real and s purious matches distributions. Because the spatial distribution of the X-ray sources in not even, we were not able to use the 
method of lDamiani et al.l d2006l) to evaluated spurious matches (which is based on the hypothesis of catalogs with an almost uniform 
spatial distribution). To see how the number of casual identifications increases with the matching radius, we produced 40 random 
X-ray catalogs by rigid translations of the X-ray coordinates (from 1' to 4' in each direction). Then, we matched with different A the 
optical-infrared catalog with these random X-ray catalogs, in order to find the distribution of casual coincidences at various A. We 
set A = 1 because this is the highest value at which the spurious matches are negligible with respect to the real ones. For an off-axis 
smaller than 3' we set the matching radius at 0".5. The number of multiple matches in our final catalog is reasonably low (860, with 
88 produced by matching the X-ray catalog). 
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